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Abstract. For any two-sided jumping a-stable process, where a e (1,2), 
we find an explicit identity for the law of the first hitting time of the origin. 
This complements existing work in the symmetric case and the spectrally one- 



sided case; cf. [32|, [15| and [28], respectively. We appeal to the Lamperti-Kiu 
representation of Chaumont et al. ll[ for real-valued self-similar Markov 
processes. Our main result follows by considering a vector-valued functional 



equation for the Mellin transform of the integrated exponential Markov ad- 
ditive process in the Lamperti-Kiu representation. We conclude our present- 
ation with some applications. 
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<N 

1 Introduction 

A Levy process is a stochastic process issued from the origin with stationary and in- 
dependent increments and cadlag paths. If X := (X t )t^o is a one-dimensional Levy 
process with law P then the classical Levy-Khintchine formula states that for all ^el, 
the characteristic exponent ^(0) := — logE(e l6IXl ) satisfies 



#(0) = ia9 + \cr 2 9 2 + f (1 - e Wx + i0a;l {W<1} )n(da;) 
2 Jr 
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where a e R, a ^ and II is a measure (the Lev?/ measure) concentrated on R\{0} such 
that L(l a x 2 )Tl(dx) < oo. 

The process (X, P) is said to be a (strictly) a-stable process if it is a Levy process 
which also satisfies the scaling property: under P, for every c > 0, the process (cX tc - a ) t ^ 
has the same law as X. It is known that a e (0,2], and the case a = 2 corresponds 
to Brownian motion, which we exclude. In fact, we will assume a e (1,2). The Levy- 
Khintchine representation of such a process is as follows: o = 0, and II is absolutely 
continuous with density given by 

c + x-( a+1 h {x>0} + c_\x\-( a+1 h {x<0} , xeR, 

where c+, c_ ^ 0, and a = (c + — c_)/(a — 1). 
The process X has the characteristic exponent 

tf(0) = c |#| Q (l-i/3tan^sgn#), 9 e R, (1.1) 

where /3 = (c + — c_)/(c+ + c_) and c = — (c+ + c_)T(— a) cos(7ra;/2). For more details, 
see Sato [30, §14]. 

For consistency with the literature that we shall appeal to in this article, we shall 
always parameterise our a-stable process such that 

sin(^ap) sin(Trap) 
c+ = 1 (a + 1) and c_ = 1 (a + 1) , 

7T 7T 

where p = P(X t ^ 0) is the positivity parameter, and p = 1 — p. In that case, the 
constant c simplifies to just c = cos(7ra(p — 1/2)). 

We take the point of view that the class of stable processes, with this normalisation, 
is parameterised by a and p; the reader will note that all the quantities above can be 
written in terms of these parameters. We shall restrict ourselves a little further within 
this class by excluding the possibility of having only one-sided jumps. This gives us the 
following set of admissible parameters: 

At = {(a, p) ■ a e (1, 2), p g (1 - 1/a, l/a)}. 

We are interested in computing the distribution of 

T = inf{t ^ : X t = 0}, 

the first hitting time of zero for X. When a > 1, this random variable is a.s. finite, while 
when a < 1, points are polar for the stable process, so T = oo a.s.; this explains our 
exclusion of such processes. This paper consists of two parts: the first deals with the case 
where X is symmetric. Here, we may identify a positive, self-similar Markov process R 
and make use of the Lamperti transform to compute T . The second part concerns the 
general case where X may be asymmetric. Here we present instead a method making 
use of the generalised Lamperti transform and Markov additive processes. 

It should be noted that the symmetric case can be deduced from the general case, 
and so in principle we need not go into details when X is symmetric; however, this case 
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provides familiar ground on which to rehearse the arguments which appear in a more 
complicated situation in the general case. Let us also note here that, in the symmetric 
case, the distribution of T has been characterised in Yano et al. [32, Theorem 5.3], and 



the Mellin transform appears in Cordero 15, equation (1.36)]; however, these authors 
proceed via a different method. 

For the spectrally one-sided case, which our range of parameters omits, the law of Tq 



has been computed by Peskir [28[. This justifies our exclusion of the one-sided case. 
Nonetheless, as we explain in Remark 13. 9[ our methodology can also be used in this 
case. 

We now give a short outline of the coming material. 

In section 121 we suppose that the stable process X is symmetric, that is p = 1/2, and 
we define a process R by 

R t = \Xt\l{t<2b}, t ^ 0, 

the radial part of X. The process R satisfies the a-scaling property, and indeed is a 
positive, self-similar Markov process, whose Lamperti representation, say £, is a Levy 
process; see section [2] for definitions. It is then known that T has the same distribution 
as the random variable 

I(ag) := exp(a^t) dt, 
Jo 

the so-called exponential functional of a£. 

In order to find the distribution of T , we compute the Mellin transform, E[/(ai£) s_1 ], 
for a suitable range of s. The result is given in Proposition 12.31 This then characterises 
the distribution, and the transform here can be analytically inverted. 

In section [3] we consider the general case, where X may not be symmetric, and our 
reasoning is along very similar lines. The process R still satisfies the scaling property, 
but it is no longer a Markov process, and the method above breaks down. However, due 



to the recent work of Chaumont et al. [11] , there is still a type of Lamperti representation 
for X, not in terms of a Levy process, but in terms of a so-called Markov additive process, 
say £. Again, the distribution of T is equal to that of I(a£) (but now with the role 
of £ taken by the Markov additive process), and we develop techniques to compute 
a vector-valued Mellin transform for the exponential function of this Markov additive 
process. Further, we invert the Mellin transform of I(cx£) in order to deduce explicit 
series representations for the law of To- 

In certain scenarios the distribution T is a very convenient quantity to have, and we 
consider some applications in section HJ for example, we give an alternative description 
of the stable process conditioned to avoid zero, and we give some identities in law similar 
to the result of Bertoin and Yor [i[ for the entrance law of a pssMp started at zero. 



2 The symmetric case 

In this section, we give a brief derivation of the Mellin transform of To for a symmetric 
stable process. As we have said, we do this by considering the Lamperti transform of the 
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radial part of the process; let us therefore recall some relevant definitions and results. 

A positive self-similar Markov process (pssMp) with self- similarity index a > is a 
standard Markov process R = (Rt)t^o with associated filtration {^t) t >o anc ^ probability 
laws (P x ) x> o, on [0, go), which has as an absorbing state and which satisfies the scaling 
property, that for every x, c > 0, 

the law of (ci? te - a ) t>0 under P x is P cx . 

Here, we mean "standard" in the sense of jij, which is to say, (<^t) t >o 1S a complete, 
right-continuous filtration, and R has cadlag paths and is strong Markov and quasi-left- 
continuous. 



In the seminal paper [23], Lamperti describes a one-to-one correspondence between 
pssMps and Levy processes, which we now outline. It may be worth noting that we have 
pre sented a slightly different definition of pssMp from Lamperti; for the connection, see 

m §o]- 

Let S(t) = ^ (i? u ) _a du. This process is continuous and strictly increasing until R 
reaches zero. Let (T(s)) sS , be its inverse, and define 

f] s = \ogR Tis) 0. 

Then rj := (f]s)s^o is a Levy process started at \ogx, possibly killed at an independent 
exponential time; the law of the Levy process and the rate of killing do not depend on 
the value of x. The real-valued process rj with probability laws (P^gr is called the Levy 
process associated to R, or the Lamperti transform of R. 

An equivalent definition of S and T, in terms of i] instead of R, is given by taking 
T(s) = ^exp^aiju) du and S as its inverse. Then, 

R t = exp(?75 W ) 

for all t ^ 0, and this shows that the Lamperti transform is a bijection. 

Let X be the symmetric stable process, that is, the process defined in the introduction 
with p = 1/2. We note that this process is such that E(e l9Xl ) = exp(— 16*1"), and remark 
briefly that it has Levy measure Il(dx) = A;|x| _a_1 dx, where 

fc = r(a + l)^. 

7T 

Connected to X is an important example where the Lamperti transform can be com- 
puted explicitly. In Caballero and Chaumont p], the authors compute Lamperti trans- 
forms of killed and conditioned stable processes; the simplest of their examples is as 
follows. Let 

t = inf{t > : X t < 0}, 

and denote by £* the Lamperti transform of the pssMp (^X t t^ t<T -^j . Then has 
Levy density 

ke x \e x - l\- [a+1 \ jel, 
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and is killed at rate k/a. 

Using this, we can analyse a pssMp which will give us the information we seek. Let 

R t = |X t |l{ t<To }, t 5= 0, 

be the radial part of the symmetric stable process. It is simple to see that this is a 
pssMp. Let us denote its Lamperti transform by £. 

In Caballero et al. [7|, the authors study the Levy process £; they find its characteristic 
function, and Wiener-Hopf factorisation, when a < 1, as well as a decomposition into 
two simpler processes when a < 1. We will now demonstrate that their expression for 
the characteristic function is also valid when a > 1, by showing that their decomposition 
has meaning in terms of the Lamperti transform. 

Proposition 2.1. The Levy process £ is the sum of two independent Levy processes, £ L 
and £ c , such that 

(i) The Levy process £ L has characteristic exponent 

^(6) -k/a, 9eR, 

where ty* is the characteristic exponent of the process £*, which is the Lamperti 
transform of the stable process killed upon first passage below zero. That is, £ L is 
formed by removing the independent killing from £*. 

(ii) The process £ c is a compound Poisson process whose jumps occur at rate k/a, 
whose Levy density is 

e y 

% c (y) = k- — -, yeR. 

Proof. Precisely the same argument as in [22[ gives the decomposition into £ L and £ c , 
and the process £ L is exactly as in that case. The particular form of the compound 
Poisson process is new. In Proposition 1 of Q, the authors decompose the Levy measure 
of £ into that of £ L and the measure with density 7r c ; they do so under the assumption 
a e (0, 1], however, this is not used in the proof, which works for any a e (0, 2). Since 
the law of a compound Poisson process is determined entirely by its Levy measure, this 
is sufficient. □ 

We now compute the characteristic exponent of £. As we have mentioned, when a < 1, 
this has already been computed in 0, Theorem 7], but whereas in that paper the authors 
were concerned with computing the Wiener-Hopf factorisation, and the characteristic 
function was extracted as a consequence of this, here we provide a proof directly from 
the above decomposition. 

Theorem 2.2 (Characteristic exponent). The characteristic exponent of the Levy pro- 
cess £ is given by 

m = 2 *n*ft-m r(V2+io/2) eeR 

KJ T(-i0/2) r((l-a)/2 + i0/2)' v ; 
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Proof. We consider separately the two Levy processes in Proposition 12. 11 
Recall that 



# L (0) 



r(a - i0)r(l + \9) 



Via) 



r(a/2 - \6)T{1 - a/2 + i0) T(a/2)r(l - a/2) 



and via the beta integral, 
^ c (0) = k 



7r (y)dy 



T(a+ 1) 



T(a/2)r(l - a/2) 



1 r(l + i0)r(a - i£ 
a 



T a + 1 



Summing these, then using product-sum identities and 8.334.2-3], 

sin(7r(a/2 — i0)J sin("7ra/2) 



#(0) = r(a-i0)r(l + i 
2 



7T 



7T 



^/ „s 7r(a — i0) — i07T 
-r(a - i0)r(l + i0) cos - sin — — 

7T 2 2 

r(a-i0) r(i + i0) 

27T- 



r(l/2 + (a - i0)/2) T(l/2 + (1 + i0)/2) r(-i0/2)r((l - a + i0)/2) ' 

Now, applying the Legendre duplication formula [17, 8.335.1] for the gamma function, 
we obtain the expression in the theorem. □ 

We now characterise the law of the exponential functional 



of the process a£, which we do via the Mellin transform, 

M(s) = K[I(aO s - 1 l 

for s e C whose real part lies in some open interval which we will specify. 

To begin with, we observe that the Laplace exponent if} of the process — a£, that is, 
the function such that Ee -20 ^ 1 = e^ z \ is given by 



. r(l/2 - az/2) r(a(l + z)/2) 
T(l/2-a(l + z)/2) T(az/2) 



Re z 6 (— 1, 1/a). 



We will now proceed via the 'verification result' [2ll . Proposition 2]: essentially, this 
result says that we must find a candidate for A4 which satisfies the functional equation 



M(s + 1) 



M(s), 



(2.2) 



for certain s, together with some additional conditions. Let us now state our result. 
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Proposition 2.3. The Mellin transform of To satisfies 

sin 7r(s-l + i) s; 



(2.3) 



for Res e (-1,2 - 1). 



Proof. Denote the right-hand side of (2.3) by f(s). We begin by noting that — a£ satisfies 
the Cramer condition tp(l/a — 1) = 0; the verification result 2l|, Proposition 2] therefore 
allows us to prove the proposition for Re s e (0, 2 — 1/a) once we verify some conditions 
on this domain. 

There are three conditions to check. For the first, we require f(s) is analytic and 
zero-free in the strip Res e (0, 2 — this is straightforward. The second requires us 



to verify that / satisfies (2.2) To this end, we expand cos and sin in gamma functions 
(via reflection formulas [17, 8.334.2-3]) and apply the Legendre duplication formula 17J, 
8.335.1] twice. 

Finally, there is an asymptotic property which is needed. More precisely, we need to 
investigate the asymptotics of /(s) as Im(s) — » go. To do this, we will use the following 
formula 



lim 



\T(x + iy)\ 



\y\ 



i 



2,71. 



(2.4) 



This can be derived from the Stirling's asymptotic formula: 



io g (r(*)) 



2 ) M 2 ) 



2: + -log(27r) + 0(2:- 1 ), 



(2.5) 



as z —*■ oo and | arg(z)| < n — 5, for fixed 5 > 0. 

Since Stirling's asymptotic formula is uniform in any sector | arg(z) | < tt — 5, it is easy 
to see that the convergence in (2.4) is also uniform in x belonging to a compact subset 
of R. Using formula (2.4) we check that |l//(s)| = 0(exp(-7r|Im(s)|)) as Im(s) — > oo, 
uniformly in the strip Re s 6 (0, 2 — 1/a). This is the asymptotic result that we require. 

The conditions of [HI, Propositon 2] are therefore satisfied, and it follows that the 
formula in the proposition holds for Re s e (0, 2 — 1/a). Since f(s) is analytic in the wider 
strip Re(s) e (—1/a, 2— 1/a), we conclude that E [/(a^) s_1 ] can be analytically extended 



to this wider strip and that identity (2.3) holds for all s e C with Re(s) e (—1/a, 2 — 1/a 
Theorem 2]). 



sec 



□ 



We note that the expression of Cordero [15j, equation (1.36)] can be derived from this 
result via the Legendre duplication formula for the gamma function. 

Now, this Mellin transform completely characterises the law of T , and we could at 
this point invert the Mellin transform to find a series expansion for the density of Tq. 
However, as we will shortly perform precisely this calculation in section [3] for the general 
case, we shall leave the Mellin transform as is and proceed to consider what happens 
when X may not be symmetric. 
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3 The asymmetric case 



With the symmetric case as our model, we will now tackle the general case where X 
may be asymmetric. The ideas here are much the same, but the possibility of asym- 
metry means that we require more complicated objects: our positive self-similar Markov 
processes become real self-similar Markov processes; our Levy processes become Markov 
additive processes; and our functional equation for the Mellin transform (12 .2p becomes 
vector-valued. 

The section is laid out as follows. We devote the first two subsections to a discussion 
of Markov additive processes and their exponential functionals, and then discuss real 
self-similar Markov processes and the generalised Lamperti representation. Finally, in 
the last subsection, we apply the theory which we have developed to the problem of 
determining the law of T for a general two-sided jumping stable process with as (1,2). 



3.1 Markov additive processes 

Let E be a finite state space and {&t)t^o a standard filtration. A cadlag process (£, J) in 
K. x E with law P is called a Markov additive process (MAP) with respect to {&t)t^o if 
( J(t)) t ^o is a continuous-time Markov chain in E, and the following property is satisfied, 
for any i e E, s, t ^ 0: 

given {J(t) = i}, the pair + s) — J{t + s)) is independent of 

and has the same distribution as — £(0), J(s)) given { J(0) = i). (3.1) 

Aspects of the theory of Markov additive processes are covered in a number of texts, 
among them [l[ and p|. We will mainly use the notation of 19], which principally works 



under the assumption that £ is spectrally negative; the results which we quote are valid 
without this hypothesis, however. 

Let us introduce some notation. We write P, = P(-|£(0) = 0, J(0) = i); and if /i is a 
probability distribution on E, we write P M = P(-|f (0) = 0, J(0) ~ //) = £ fcJi . /x(«)P*. We 
adopt a similar convention for expectations. 



It is well-known that a Markov additive process (£, J) also satisfies (3.1) with t replaced 



by a stopping time. Furthermore, it has the structure given by the following proposition; 



see [2|, §XI.2a] and [l9j, Proposition 2.5]. 



Proposition 3.1. The pair (£, J) is a Markov additive process if and only if, for each 
i,j e E, there exist a sequence of iid Levy processes (Ci)n^o Q,nd a sequence of iid random 
variables (U-j) n ^o, independent of the chain J, such that if T = and (T n ) n5 >i are the 
jump times of J, the process £ has the representation 

£(t) = l {n>0} (£(T n -) + U] {Tn _ u(Tn) ) + Cj {T Jt - T n ), t e [T n , T n+1 ), n>0. 

For each i 6 E, it will be convenient to define, on the same probability space, as 
a Levy process whose distribution is the common law of the processes in the above 
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representation; and similarly, for each i,j e E, define to be a random variable having 
the common law of the UTi variables. 

Let us now fix the following setup. Firstly, we confine ourselves to irreducible Markov 
chains J. Let the state space E be the finite set {1, . . . , N}, for some iV e N. Denote 
the transition rate matrix of the chain J by Q = (qij)ijeE- F° r each i e E, the Laplace 
exponent of the Levy process & will be written tpi, in the sense that e^ z > = E(e*^ i ' 1 ^), 
for all z e C for which the right-hand side exists. For each pair of i,j e E, define the 
Laplace transform Gij(z) = ~K(e zUl1 ) of the jump distribution U^, where this exists; write 
G(z) for the N x N matrix whose (i, j)th element is Cry (z) .We will adopt the convention 
that Uij = if Qij = 0, i j, and also set Uu = for each i e E. 

A multidimensional analogue of the Laplace exponent of a Levy process is provided 
by the matrix-valued function 



for all z e C where the elements on the right are defined, where o indicates elementwise 
multiplication, also called Hadamard multiplication. It is then known that 



for all z e C where one side of the equality is defined. For this reason, F is called the 
matrix exponent of the MAP £. 

We now describe the existence of the leading eigenvalue of the matrix F, which will 
play a key role in our analysis of MAPs. This is sometimes also called the Perron- 
Frobenius eigenvalue; see @, §XI.2c] and [19, Proposition 2.12]. 

Proposition 3.2. Suppose that z e C is such that F(z) is defined. Then, the matrix 
F(z) has a real simple eigenvalue k(z), which is larger than the real part of all its other 
eigenvalues. Furthermore, the corresponding right-eigenvector v(z) may be chosen so 
that Vi(z) > for every i e E, and normalised such that 



where n is the equilibrium distribution of the chain J. 

This leading eigenvalue features in the following probabilistic result, which identifies a 
martingale (also known as the Wald martingale) and change of measure analogous to the 
exponential martingale and Esscher transformation of a Levy process; cf. 0, Proposition 
XI.2.4, Theorem XIII.8.1]. 

Proposition 3.3. Let 



F(z) = diag(^), . . . , ij N (z)) + Q o G(z) 



(3.2) 



E i (e z ^;J(t)=j) = (e F ^ t ) 



h j e E, 



TTv(z) = 1 



(3.3) 



M(t, 7 ) =e 



wj(o)(7) 



t > 0, 



for some 7 such that the right-hand side is defined. Then, 
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(i) M(-, 7 ) is a unit-mean martingale with respect to (@t)t>o under any initial distri- 
bution o/(f(0),J(0)). 

(ii) Define the change of measure 



= M(t, 7 ). 

Under PW, process £ is still a Markov additive process, and it has the following 
characteristics, for each i,j e E: 

• PM(*7« e dx) = 7r ^ T P(f/ ii e dx), and Zience G^(«) = ^ . , , 

• $ } = ^j^ G M *nd 

Furthermore, 

= (diag(^( 7 ), i e E))- l [F{z + 7) - «( 7 )Id] diagfafr), 2 e E), 

and hence, 

k^\z) = k(z + 7) — k( 7 ). 

Making use of this, the following proposition can be obtained; the properties of k 
given here are often used in the literature, but for convenience, we also provide a short 
proof. 

Proposition 3.4. Suppose that F is defined in some open interval D o/R. Then, the 
leading eigenvalue k of F is smooth and convex on D. 



Proof. Smoothness follows from results on the perturbation of eigenvalues; see [19|, Pro- 
position 2.13] for a full proof. 

To prove that k is convex, we begin by proving that for A 6 (0, 1) such that z, Xz e D, 

k(Xz) Xk(z). (3.4) 



Taking expectations in the martingale from Proposition 13. 3[ and using the normalisation 



condition (3.3) , we see that, starting J in equilibrium, 

E n [e z ^v J{1) (z)] = e^\ 

Now let p = 1/X > 1, and let q be its conjugate exponent, that is, - + - = 1. We apply 
Holder's inequality as follows. 

e «(./p) =E K [eW>M% JW {z)] < (Eje^( 1 )^ ) ( 2 )]) 1/p (E.[^ 1) ( 2 )]) 1 ^ 

= {E n [e^% J(1) (z)]) 1/p (nv(z))^ 
= (E n [e^% J(1) (z)]) 1/p = e «W/" 
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Replacing 1/p by A, we obtain (3.4) 



Now let u,v e D and A 6 (0,1); recall from Proposition 13.31 that for the leading 
eigenvalue under the change of measure, k^ u '(z) = k(z + u) — k(u). We then calculate: 

«((1 - X)u + Xv) - k{u) = k (u) (X(v - u)) 

< Xk {u) (v -u) 
= X(k(v) — n(u)) 

and so 

k((1 - X)u + Xv) < (1 - X)k(u) + Xk(v), 
which completes the proof. □ 

3.2 The Mellin transform of the exponential functional 

In section [2J we studied the exponential functional of a certain Levy process associated 
to the radial part of the stable process; now we are interested in obtaining some res- 
ults which will assist us in computing the law of an integrated exponential functional 
associated to Markov additive processes. 
For a MAP f , let 

/(-$= exp(-e(t)) df. 
Jo 

One way to characterise the law of is via its Mellin transform, which we write as 

Ai(s). This is the vector in ~R N whose ith element is given by 

Mi(s) =E i [I(-£) s - 1 ], ieE. 

We will shortly obtain a functional equation for Ai, analogous to the functional equa- 



tion (2.2) which we saw in section [2j For Levy processes, proofs of the result can be 
found in Carmona et al. 0, Proposition 3.1], Maulik and Zwart 25, Lemma 2.1] and 
Rivero 0, Lemma 2]; our proof follows the latter, making changes to account for the 
Markov additive property. 

We make the following assumption, which is analogous to the Cramer condition for 
a Levy process; recall that k is the leading eigenvalue of the matrix F, as discussed in 
section 13. 11 

Assumption 3.5 (Cramer condition for a Markov additive process). There exists zq < 
such that F(s) exists on (zq,0), and some 9 e (0, — zq), called the Cramer number, such 
that k(-0) = 0. 

Since the leading eigenvalue k is smooth and convex where it is defined, it follows 
also that k{— s) < for s e (0, 9). In particular, this renders the matrix F(—s) negative 
definite, and hence invertible. Furthermore, it follows that re'(0— ) > 0, and hence (see 
0, Corollary XI. 2. 7] and 0, Lemma 2.14]) that £ drifts to +00 independently of its 
initial state. This implies that J(— £) is an a.s. finite random variable. 
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Proposition 3.6. Suppose that £ satisfies the Cramer condition (Assumption 1 3. 5\) with 
Cramer number 9 e (0, 1). Then, M.(s) is finite and analytic when Re s e (0, 1 + 9), and 
we have the following vector-valued functional equation: 

M(s + l) = -s{F{-s))- l M{s), se(0,9). 

Proof. At the end of the proof, we shall require the existence of certain moments of the 
random variable 



Qt 



du, 



and so we shall begin by establishing this. 

Suppose that s e (0, 9], and let p > 1. Then, by the Cramer condition, it follows that 
k(—s/p) < 0, and hence for any u ^ 0, e~ UK (~ S//p ) ^ 1. 

Recall that the process 

M(u,z) = e ^H-«W«^44. u >® 

Vj(0){Z) 

is a martingale (the Wald martingale) under any initial distribution (£(0), J(0)), and set 

Ffz) = minvJz) > 0, 

so that for each j e E, Vj(z)/V(z) ^ 1. 

We now have everything in place to make the following calculation, which uses the 
Cramer condition in the third line, and the Doob maximal inequality in connection with 
the Wald martingale in the fourth line. 



Ei[Ql] < fE t 



SMv[M(u,-s/p) Vl (-s/p)(V(-s/p))- l Y 



< t s Vi{-s/p) p v{~s/pY p ( — ^— j- j Ei[M(t,-s/ P y] 



p 



p — 1 



p 



-tpK( — s/p) 



< 00. 



Now, it is simple to show that for all s > 0, t ^ 0, 



(f 



£(«)( e -(€(«+«)-£(u)) dv 



dw. 



For each i e E, we take expectations and apply the Markov additive property. 



-^du) 



E j 



2 f E i \e- s ^;J(u)=j 



dt; 



d-u 
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Since < s < 9 < 1, it follows that 1 1 x | s — \y\ s \ ^ \x — y\ s for any x, y e R, and so we see 
that for each i e E, the left-hand side of the above equation is bounded by Ki(Qf) < go. 
Since (e F ^ u ) .. ^ 0, it follows that E i [/(-^) s - 1 ] < oo also. 

If we now take t —*■ oo, the left-hand side is monotone increasing, while on the right, the 
Cramer condition ensures that F(—s) is negative definite, which is a sufficient condition 
for convergence, giving the limit: 

M{s + 1) = s e (0, 9). 

Furthermore, as we know the right-hand side is finite, this functional equation allows 
us to conclude that Ai(s) < oo for all s e (0, 1 + 9). It then follows from the general 
properties of Mellin transforms that Ai(s) is finite and analytic for all s s C such that 
Rese (0,1 + 9). 

□ 



3.3 Real self-similar Markov processes 

In section [5J we studied a Levy process which was associated through the Lamperti 
representation to a positive, self-similar Markov process. Here we see that Markov 
additive processes also admit an interpretation as Lamperti-type representations of real 
self-similar Markov processes. 

The structure of real self-similar Markov processes has been investigated by Chy- 
biryakov 14] in the symmetric case, and Chaumont et al. [11] in general. Here, we 



give an interpretation of these authors' results in terms of a two-state Markov additive 
process. We begin with some relevant definitions, and state some of the results of these 
authors. 

A real self-similar Markov process with self- similarity index a > is a standard (in 
the sense of 0]) Markov process X = (X t ) t>0 with probability laws (P x )xeM\{o} which 
satisfies the scaling property, that for all x e M\{0} and c > 0, 

the law of (cX tc - a ) t>0 under is P cx . 



In [11] the authors confine their attention to processes in 'class C.4'. An rssMp X is 
in this class if, for all x / 0, P x (3t > : X t X t - < 0) = 1; that is, with probability one, 
the process X changes sign infinitely often. As with the stable process, define 

T = inf{t ^ : X t = 0}. 

For such processes, X may be identified up to the time T as the time-changed expo- 
nential of a certain complex- valued process £, which we call the Lamperti-Kiu repres- 



entation of X. The following result is taken from [11 . 



Proposition 3.7 ([11, Theorem 6]). Let X be a rssMp in class C-4, and let x 0. It 

is possible to define independent sequences (£ ±,fc )fc>o> (C^Wcb {U ±,k )k^o of iid random 
objects with the following proprties: 
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1. The elements of these sequences are distributed such that: the ^- are real-valued 
Levy processes; are exponential random variables with parameters q ± ; and U ± 
are real-valued random variables. 



2. For each i^O, define the following objects 
fcW C W Z7(*)) = 



(e + ' fc ,C + ' fc ,^ + ' fc ), if sgn(x)(-l) fc = 1 
{£-' h ,(-> k ,U-> k ), if sgn(x)(-l) fc = -1, 



n-1 

r = o, r n = Y J C (k \ 

N t = max{n ^ : % < *}, 
o~t = t — 7jvt> 

ATt-l 

^(fc) 

fc=0 



£ = ^ + 2 (45, + ^ (fc) ) + wrJVi, t 2* 0, 

fc=0 

r(t) = inf|s > : J |exp(a£ M )| du > t\ 



M^j, f<T . 
Then, the process X under the measure P x has the representation 

X t = xexp(£ T{t )), 0^t<T . 

The abundance of notation necessary to be precise in this context may obscure the 
fundamental idea, which is as follows. At any given time, the process £ evolves as a 
Levy process moving along a line Imz = nN, up until an exponential 'clock' 
(corresponding to the process X changing sign) rings. At this point the imaginary part 
of £ is incremented by %, the real part jumps by U ± , and the process begins to evolve 
as the other Levy process, £ + . 

Particularly in light of the above discussion, our interpretation is simple to state. 

Proposition 3.8. Let X be an rssMp, with Lamperti-Kiu representation £. Define 
furthermore 

1, if n is odd 

H = { 

12, if n is even. 

Then, for each i^O, the process 

(£(*), J(t)) = (Re£ t , [lm£ t /ir + l {x >o}]) 
defined in Proposition \3. 7| is a Markov additive process with state space E = {1,2} ; and 



X under Y x has the representation 

X t = zexp(£(r(*)) + in(J(r(t)) + 1)), < t < T , 
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where we note that (£(0), J(0)) is equal to (0, 1) if x > 0, or (0, 2) if x < 0. Furthermore, 
the time- change r has the representation 



r(t) = inf < s > : exp(a£(u)) du > £|x| 



t<T , 



(3.5) 



m terms of the real-valued process £. 



Note that the MAP (£, J) under Pi corresponds to the rssMp X started at a point 
x > 0, and the MAP under P 2 corresponds to the rssMp started at a point x < 0. 
Furthermore, we observe from the form (3.5) of the time-change r that under P x , for 



any x 0, the following identity holds for T , the hitting time of zero: 



\x\ 



f 

Jo 



; Q?(u) du. 



Implicit in this statement is that the MAP on the right-hand side has law Pi if x > 0, 
and law P 2 if 2 < 0. This observation will be exploited in the coming section, in which 
we put together the theory we have outlined so far. 



3.4 The hitting time of zero 

We now return to the central problem of this paper: computing the distribution of T 
for a stable process. We already have in hand the representation of T for an rssMp as 
the exponential functional of a MAP, as well as a functional equation for this quantity 
which will assist us in the computation. 

Let X be the stable process with parameters (a, p) e A s t, defined in the introduction. 
We will restrict our attention for now to X under the measures P+i; the results for other 
initial values can be derived via scaling. 

Since X is an rssMp, it has a representation in terms of a MAP (£, J); futhermore, 
under P+i, 



T 



f 

Jo 



3 <*M ds = I(aO; 



to be precise, under Pi the process £ is under Pi, while under P_i it is under P 2 



In 11 



^4.1], the authors calculate the characteristics of the Lamperti-Kiu represent- 
ation for X, that is, the processes £ ± , and the jump distributions U ± and rates q ± . 
Using this information, and the representation (3.2) , one sees that the MAP (— a£, J) 
has matrix exponent 



/ r(a(l + z))r(l - az) 



F(z) 



V 



T(ap + az)T(l — ap — az) 
T{a(l + z))T(l - az) 

r(ap)r(l — ap) 



r(q(l + z))r(l - az) \ 

r(ap)r(l — ap) 
r(a(l + z))T(l - az) 

T(ap + az)T(l — ap — az) ) 



for Re z e (—1, 1/a). 
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Remark 3.9. It is well-known that, when X does not have one-sided jumps, it changes 
sign infinitely often; that is, the rssMp X is in [ll|'s class C.4. When the stable process 
has only one-sided jumps, which corresponds to the parameter values p = 1 — l/a, l/a, 
then it jumps over at most once before hitting it; the rssMp is therefore in class 
C.l or C.2 according to the classification of [11]. The Markov chain component of the 
corresponding MAP then has one absorbing state, and hence is no longer irreducible. 
Although it seems plain that our calculations can be carried out in this case, we omit it 



for the sake of simplicity. As we remarked in the introduction, it is considered in [28 



We now analyse F in order to deduce the Mellin transform of Tq. By considering the 
equation det F(z) = with z = l/a — 1, it is not difficult to deduce that n(l/a — 1) = 0; 
that is, — a£ satisfies the Cramer condition (Assumption 13. 5j) with Cramer number 9 = 
1 - l/a. 
Define 

Ms) := EipT 1 ] = Eif/Kr 1 ], Ma) := E_i[T s_1 ] = E^K) 5 " 1 ], 

which by Proposition 13.61 are defined when Res e (0,2 — l/a). This proposition also 
implies that 



B( S ) 



" fi(s + l) ' 




" fi(s) ' 


_ M* + 1) _ 




_ Ms) _ 



s 6 (0,1 - l/a), 



(3.6) 



where B(s) := —F(—s)/s. Using the reflection formula for the gamma function we find 
that 



a 

B(s) = — T(a — as)T(as) 

TV 



sin(7ra(/3 — s)) — sin(7rap) 
— sin(7rap) sin(7roi(p — s)) 



Res e (—l/a, 1), s^O, 



Using the above formula and Ptolemy's identity 

sin(5i + 5 2 ) sin(<5 2 + S3) = sin(5i) sin(5 3 ) + sin(<5i + 5 2 + S 3 ) sm(5 2 ) 
with 81 = irap, 5 2 = —nas and £3 = nap, we conclude that 

Via — as)T(as) 



det(B(s)) = —a „ 

y y " r(l - a + as)Y(l - as) 

Therefore, if we define A(s) = (B(s)) -1 , we have 



Re s g (-l/a, 1), s ^ 0. 



(3.7) 



{3A 



sin(7ra(p — s)) sin(7rap) 
sin(7rap) sin(7ra(p — s)) 



A(s) = r(l - a + as)lYl - as) 

ixa 

for Re s e (1 — 2/a, 1 — l/a), and may rewrite (3.6) in the form 

= A(s)x , se (0,1 -l/a). 

_ J2{S) 

The following theorem is our main result. 



Ms + i) 



(3.9) 
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Theorem 3.10. For —1/a < Re(s) < 2 — 1/a we have 



F IT*' 1 ! Sin ^ Sin - a + ^ r (! + « - Qg ) mm 

U J sin(Trp) sin(*(l-a + as)) T(2 s) ' { } 



The corresponding expression for E_i[T s x ] can be obtained from (3.10) by changing 
p^p. 



Let us denote the function in the right-hand side of (3.10) by h\(s), and by /^(s) the 
function obtained from hi(s) by replacing p p. Before we are able to prove Theorem 
13.101 we need to establish several properties of these functions. 

Lemma 3.11. 

(i) There exists e e (0, 1 — 1/a) such that the functions hi(s), h 2 (s) are analytic and 
zero-free in the vertical strip < Re(s) < 1 + e. 

(ii) For any —cc<a<b<+cc there exists C > such that 

e -7r|Im«l < |^( s )| < e -%(a-l)\1m(s)\ ^ { = ^ 3 

for all s in the vertical strip a < Re(s) < b satisfying |Im(s)| > C. 

Proof. It is clear from the definition of h\{s) that it is a meromorphic function. Its 
zeroes are contained in the set 

{2, 3, 4, . . . } u {1 - 1/a + n/(ap) : n e Z, n ^ 0} 

and its poles are contained in the set 

{1 + n/a: n ^ 1} u {n — 1/a : neZ,n^l}. 

In particular, there are no zeroes or poles of hi(s) in the strip < Re(s) < 1. The same 
is clearly true for h 2 (s), which proves part (i). 



We now make use of Stirling's formula, as we did in section [2J Applying (2.4) to hi(s) 
we find that as s —*■ oo (uniformly in the strip a < Re(s) < b) we have 

log(|/n(s)|) = -|(1 + a 2ap)|Im(s)| + 0(log(|Im(a)|)). 

Since for a > 1, the admissible parameters a, p must satisfy a — 1 < ap < 1, this shows 
that 

a — 1 < 1 + a — lap < 3 — a < 2, 
and completes the proof of part (ii). □ 



Lemma 3.12. The functions h\{s), h%{s) satisfy the system of equations (3.9) 
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Proof. Denote the elements of the matrix A(s) by Ay(s). Multiplying the first row 
of A(s) by the column vector [h^s) , h 2 (s)] T , and using identity sin(7rp) = sin(7rp), we 
obtain 



A ll (s)h 1 (s) + A l2 (s)h 2 (s) 

1 sin (=) T(l - as) 



7ra sin(7T/5) sin (^(1 — a + as)) 



r(l + a 



as} 



r(2 - s) 



T(l — a + as) 



x | sin(7ra(p — s)) sin (7rp(l — a + as)) + sin(7rap) sin (np(l — a + as)) }. 

Applying identity T(z + 1) = zY(z) and reflection formula for the gamma function, we 
rewrite the expression in the square brackets as follows: 



T(l + a — as) 

r(2 - s) 



r(l — a + as) 



ar(a — as) 

r(i - *) 



T(l — a + as) 



7ra 



sin(7ra(l — s))T(l — s) 



Using Ptolemy's identity (3.7) with 5i = ira(l — s), 8 2 = 7ra(s — p) and £3 = 7r(p — as + 
cups), we obtain 

sin(7ra(p — s)) sin (7rp(l — a + as)) + sin(7rap) sin (np(l — a + as)) 

= sin(7ra(l — s)) sin(7rp(l + as)). 

Combining the above three formulas we conclude 

sin (^) sin (np(l + as)) T(l - as) 



A 11 (s)/i 1 (s) + A 12 (s)h 2 (s) 



h(s + l). 



sin(7rp) sin (^(1 - a + as)) T(l - s) 

The derivation of the identity A 2 i(s)hi(s) + A 22 (s)h 2 (s) = h 2 (s + 1) is identical. We 
have now established that two functions hi(s) satisfy the system of equations (3.9) □ 

Proof of Theorem \3.1(K Our goal now is to establish the uniqueness of solutions to the 
system (3.9) in a certain class of meromorphic functions, which contains both hi(s) and 
fi(s). This will imply /ij(s) = fi(s). Our argument is similar in spirit to the proof of 



Proposition 2 in [21 



First of all, we check that there exists e e (0, 1/2 — l/(2a)), such that the functions 
/i(s), /2(s) are analytic and bounded in the open strip 

S e = {s e C : e < Re(s) < 1 + 2s} 



This follows from Proposition 13.61 and the estimate 

\fi(s)\ = IEiPT 1 ]! < E^TT 1 !] = E^T*^- 1 ] = A(Re(s)). 

The same applies to f 2 . Given results of Lemma I3.11[ we can also assume that e is small 
enough, so that the functions /ij(s) are analytic, zero- free and bounded in the strip S £ . 

Let us define D(s) := fi(s)h 2 (s) — f2(s)hi(s) for s e S E . From the above properties 
of fi(s) and /ij(s) we conclude that D(s) is analytic and bounded in S £ . Our first goal 
is to show that D(s) = 0. 
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If both s and s + 1 belong to S £ , then the function D(s) satisfies the equation 

1 r(l - a + as)T(l 



D s + 1 



as} 



-D is) 



a c 



T(a — as)T(as) 

as is easily established by taking determinants of the matrix equation 



(3.11) 



/ 2 (s + l) h 2 {s + l) 



Ms) 



f2(s) 



ht(s) 
h 2 (s) 



and using (3.8) and the identity A(s 
Define also 

„. , T(s-l)T(a 
G(s) : 



B(s 



,-i 



as) 



r(l — sjIY— a + as) 



Sill 7T 



1 

s + - 

a 



It is simple to check that: 



(i) G satisfies the functional equation (3.11 



(ii) G is analytic and zero- free in the strip S £ ; 

(iii) |G(s)| — > oo as Im(s) —*■ oo, uniformly in the strip S £ (use (2.4)). 



We now define H(s 
is analytic in the strip S E , while property 
both s and s + 1 belong to S £ . Therefore, we can extend H(s) to an entire function 
satisfying H(s + 1) 



D(s)/G(s) for s e S £ . The property (ii) guarantees that H 

H(s) if 



and (3.11) show that H(s + 1) 



m 



H(s) for all seC. Using the periodicity of H(s), property 
of the function G(s) and the fact that the function D(s) is bounded in the strip <S £ , we 
conclude that H(s) is bounded on C and H(s) — > as Im(s) —> oo. Since H is entire, it 
follows that H = 0. 

So far, we have proved that for all s e S £ we have fi(s)h,2(s) = f2(s)hi(s). Let us 
define iu(s) := fi(s)/hi(s) = f2(s)/h 2 (s). Since both fi(s) and /ij(s) satisfy the same 



functional equation (3.9) , if s and s + 1 belong to iS e we have 



+ + 1) =/i(a + l) 

= A u (a)/ 1 (s)+A 12 ( S )/ 2 ( S ) 

= w(s)[A u (s)hi(s) + A 12 (s)/i 2 (s)], 

and therefore w(s + 1) = w(s). Using again the fact that f\ and hi are analytic in this 
strip and hi is also zero free there, we conclude that w(s) is analytic in <S £ , and the 
periodicity of w implies that it may be extended to an entire periodic function. Lemma 
I3.11( ii) together with the uniform boundedness of fi(s) in S £ imply that there exists a 
constant C > such that for all s e S £ , 

\w(s)\ < Ce wllm{s)l . 

By periodicity of to, we conclude that the above bound holds for all s e C. Since w 
is periodic with period one, this bound implies that to is a constant function (see the 
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proof of Proposition 2 in [21]). Finally, we know that = hi(l) = 1, and so we 

conclude that w(s) = 1. Hence, fi(s) = hi(s) for all s s S £ . Since h,{s) are analytic in 
the wider strip < Re(s) < 2 — 1/a, by analytic continuation (see |24j . Theorem 2]) 

we conclude that (3.10) holds for all s in — 1/a < Re(s) < 2 — 1/a. 



□ 



Remark 3.13. Since the proof of Theorem 13.101 is based on a verification technique, it 



does not reveal how we derived the formula on the right-hand side of (3.10), for which 



we took a trial and error approach. As the expression in (3.10) is already known for 
spectrally one-sided case (p = 1/a or p = 1 — 1/a) due to [28l|. and for the symmetric 
case (p = 1/2) due to [32|, [15|, we sought a function which interpolated these three cases 



and satisfied the functional equation (3.9) After a candidate was found, we verified that 
this was indeed the Mellin transform, using the argument above. 



We turn our attention to computing the density of T . Let us define 
and 



mm neZ \x-n\ 



u 



{isl: lim - log 1 1 = 0}). 



This set was introduced in [18|, where it was shown that £ is a subset of the Liouville 
numbers, which implies x e £ if and only if the coefficients of the continued fraction 
representation of x grow extremely fast. It is known that C is dense, yet it is a rather 
small set: it has Hausdorff dimension zero, and therefore its Lebesgue measure is also 
zero. 

For aeiwe also define 



■ — < V( V 2)log(.Y -2))}. 



/C(a) = {N e N : || (N — f > exp 
Proposition 3.14. Assume that a ^ Q. 

(%) The set JC(a) is unbounded and has density equal to one: 

#{/C(a)n[l,n]} 



lim 

n— >oo 



1. 
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(ii) If a $ C, the set N\/C(a) is finite. 

Proof. Assume that e e (0, 1/2). It is well-known that for any irrational a the sequence 
|| (N — |)q;|| is uniformly distributed on the interval (0, 1/2). Therefore, the density of 
the set of all iV such that 



\\(N 



)a\\ > e 



is equal to 1 — 2e. At the same time, for all iV large enough, we will have exp( 
2) log(iV — 2)) < e. Thus the density of the set of all iV such that 

|| (N |)a|| > exp(-^i(iV - 2) log(iV - 2)) 



a-l 



(N- 
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is greater than 1 — 2s, which proves part (i). 

Assume that a $ C Since lim n ^ +00 ^ log||na|| = 0, there exists C > such that for 
all n we have \\na\\ > C2~ n . Then for all N we have 



\\(N - ±)a\\ > -\\(2N - l)a\\ > C2~™ . 
Since for all N large enough it is true that 

C2- 2N > exp(-^(iV - 2) log(iV - 2)), 

we conclude that all N large enough will be in the set /C(a), therefore the set N\K,(a) 
is finite. □ 

Theorem 3.15. Let p be the density of Tq under Pi. 

(i) If a ^ <Q then for all t > we /iove 

sin (-J ^— , _ , , sin 



p(t) 



lim 

NeK(a) 



7TSin(7T/3) 



/ j sm(7rp(k + 1)) 



Kfc<a(iV-i)-l 



sin (*(* + !)) 



r (- + 1 ), « 

k\ { ' 



fc -i r i- 



sin(£) ^ ™M£(HI r w+i 



7rsin(7rp) sin(7raA;) T (afc - 1) 

The above limit is uniform for t e [e, oo) and any s > 0. 



(3.12) 



(mJ If a = m/n (where m and n are coprime natural numbers) then for all t > we 
have 



p(t) 



•rn(-) sin(^) r(^ + l 

. y»{. V sinvrpfc + l . U ; r U ; , 

sin vrp ^ v PV ; sin(^ (fc + 1)) fc! 



sin 

7rsin(7rp) , t 

k^ — l (mod m) 



y sm(7capk) r (k - Q t _ k _ u 



7rsin(7rp) sin(7ra&) r (a/c — 1) 

fc^O (mod n) 

Sil1 (a) 2 V, . ' V«" -/ p ,. w -,.„-i+± 



?r 2 a sin(Trp) ^ 



(3.13) 



where 



Rk(t) := 7vapcos(7rpkm) 

— sin(7rpfcm) [7rcot (-) — ^f&n — -) + aip{km — 1) + log(t)] 



and -0 is i/ie dig amma function. The three series in (3.13) converge uniformly for 
t e [s, oo) and any s > 0. 
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(Hi) For all values of a and any c > 0, the following asymptotic expansion holds as 
t | 0: 



Proof. Recall that hi(s) = Ei[Tg x ] denotes the function in (3.10) According to Lemma 



I3.11l (ii). for any x e K, + iy) decreases to zero exponentially fast as y — » go. This 
implies that the density of To exists and is a smooth function. It also implies that p(t) 
can be written as the inverse Mellin transform, 

p(t) = — f h^r'ds. (3.14) 



2vri 

1+iR 

The function hx(s) is meromorphic, and it has poles at points 

{s$ := 1 + n/a : n ^ 1} u {s^ 2) := n - 1/a : n > 2} u {4 3) : = -w - 1/a : ra > 0} 

If a ^ Q, all these points are distinct and h 1 (s) has only simple poles. When a e Q, 
some of Sn and Sm will coincide, and hi(s) will have double poles at these points. 

Let us first consider the case a $ Q, so that all poles are simple. Let iV e K.(a) and 
define c = c(N) = N + 5 — -. Lemma f3.1 llf ii) tells us that decreases exponentially 
to zero as Im(s) — * 00, uniformly in any finite vertical strip. Therefore, we can shift the 
contour of integration in |(3.14) and obtain 



pit) = -2Res s=s ( 1) (/ il ( s )r s ) -^Res^ih^r*) 

n m 

+ J htft-'ds, (3.15) 

c(A0+iR 

where 2 n anc ^ 2m indicate summation over n ^ 1 such that < c(N) and over m ^ 2 
such that Sm < c(iV), respectively. Computing the residues we obtain the two sums in 



the right-hand side of (3.12) 



Now our goal is to show that the integral term 

2m J 



c(N)+iR 

converges to zero as iV —*■ +00, N e JC(a). We use the reflection formula for the gamma 
function and the following inequalities 

|sin(7rx)| > ||x||, |sin(x)| cosh(y) < |sin(a; + iy)\ < cosh(?/), x,yeW, 
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to estimate hi(s), s = c(N) + iu, as follows 



Ms) 



sin (tt/3(1 — a + as)) sin(7rs) T(a(s — 1)) 



sin (^(1 — a + as)) sin(7ra(s — 1)) T(s — 1) 



Ci cosli(7rapn) 



|a(iV— |)|| cosh(7rctu) 
Using Stirling's formula ( 12. 5p . we find that 



r(a(s-l)) 



T(s-l) 



(3.16) 



r(g) 



ae 



-s((a-l) log(s)+ J 4)+0(s- 1 ) 



s — ► oo, Re(s) > 0, 



where A := 1 — a + alog(a) > 0. Therefore, there exists a constant C 2 > such that 
for Re(s) > we can estimate 



T(s) 



T(as) 



< C 2 |e _s((a_1)log(s)+j4)) | 

= <7 2e -(a-l)Mslog( S ))-ARe( S )) 

£j g— (a— l)Re(s) log \s\— 4Re(s)+(a- 1) Im(s) arg(s) 

< (7 2e -( a - 1 )R e W 1 °g(M«))+("-i)|im(s)|f _ 



Combining the above estimate with (3.16) and using the fact that N e K{a) we find 
that 



\hi{c{N) + \u)\ < 



CiC 2 cosh.{napu) _ (a _ 1 



\a(N — |)|| cosh(7rcra) 



(a-l)( C (JV)-l)log(c(AO-l) + (a-l)|u|§ 



< CiC 2 e 2 



KiV-2) log(7V-2) 



coshjnapu) c(a _ 1)M 
cosh(7raw) 



Note that the function in the right-hand side of the above inequality decreases to zero 
exponentially fast as \u\ — > oo (since ap + |(a; — 1) — a < 0), and hence in particular is 
integrable on R. Thus we can estimate 



\In®\ 



f -c(N) 

2tt 



f hiciN) + iu))r h 



dit 



< 



f -c(N) 

2vr 



| /H (c(JV) +iw)) | du 



,-c(JV) 
Z7T 



a=l(JV-2)log(JV-2) 



co^(7ra^0 e(a _ 1)Mf du 
cosh(nau) 



When N — * oo, the quantity in the right-hand side of the above inequality converges to 
zero for every t > 0. This ends the proof of part (i). 

The proof of part (ii) is very similar, and we offer only a sketch. It also begins with 



(3.15) and uses the above estimate for hi(s). The only difference is that when a e Q 
some of Sn^ and Sm will coincide, and hi(s) will have double poles at these points. The 
terms with double poles give rise to the third series in (3.13) In this case all three series 
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are convergent, and we can express the limit of partial sums as a series in the usual 
sense. 

The proof of part (iii) is much simpler: we need to shift the contour of integration in 
(3. 15)| in the opposite direction (c < 0). The proof is identical to the proof of Theorem 
9 in 0. □ 



The next corollary shows that, for almost all irrational a, the expression (3.12) 
be written in a simpler form. 



can 



Corollary 3.16. //a^£uQ then 



pit) 



TrsinM^ 1 ^ "sin (£(* + !)) 



fe! 



sin (|) y sin(Trapfc) T (k — ~) ^_ fc _ 1+ _ 
7rsin(7r/)) ~ sm(irak) T (ak — 1) 



(3.17) 



The two series in the right-hand side of the above formula converge uniformly for t e 
[e, oo ) and any e > 0. 

Proof. As we have shown in Proposition ^. 14\ if a $ C u Q then the set N\/C(a) is finite. 



Therefore we can remove the restriction iV e K(ot) in (3.12) 



The remaining part of the proof if based on the following simple result: If the two 
series J] a n and J] 6 n converge, then for any increasing, unbounded sequences c n and <i n 



lim 



n=l n=l 



Therefore, to finish the proof of Corollary 13.161 it is enough to show that both series in 



(3.17} converge. 

In jlil. Proposition 1] it was shown that x e C iff x^ 1 e C. Therefore, according to 
our assumption, both a and 1/a are not in the set C. From the definition of C we see 
that there exists C > such that \\ocn\\ > C2~ n and ||a~ 1 77.|| > C2~ n for all integers n. 



Using the estimate |sin(7nc)| ^ \\x\\ and Stirling's formula (2.5), it is easy to see that 



both series in in (3.17) converge (uniformly for t e [e, oo) and any e > 0), which ends 



the proof of the corollary. 



□ 



Remark 3.17. Note that formula (3.17) may not be true if a e C, as the series may fail 



to converge. An example where this occurs is given after Theorem 5 in 20]. 



4 Applications 

4.1 Conditioning to avoid zero 



In ll|, §4.2], Chaumont, Panti and Rivero discuss a harmonic transform of a stable 
process with a > 1 which results in conditioning to avoid zero. The results quoted in 
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that paper are a special case of the notion of conditioning a Levy process to avoid zero, 
which is explored in Pantf 26|, |27 . 

In these works, in terms of the parameters used in the introduction, the authors define 



h(x) 



T(2 - a) sin(7ra/2) 



j3 sgn(a;))|a; 



a-l 



x e 



(4.1) 



C7r(a-l)(l + /3 2 tan 2 (7ra/2)) v 
If we write the function h in terms of the (a, p) parameterisation which we prefer, this 



gives 



h(x) 



x > 0, 



7T 



and the same expression with p replaced by p when x < 0. 

In 27J, Panti proves the following proposition for all Levy processes, and x e R, 
with a suitable definition of h. Here we quote only the result for stable processes and 
i^O. Hereafter, {<&t)t>o * s the standard filtration associated with X, and n refers to 
the excursion measure of the stable process away from zero, normalised (see 27J, (2.14)] 
and [16, (4.11)]) such that 



n(l 



l/li«(0), 



where ( is the excursion length and u q is the g-potential density of the stable process. 
Proposition 4.1 ( [271 . Theorem 2.2, Theorem 2.6]). Let X be a stable process, and h 



the function in (4-1) 



(i) The function h is invariant for the stable process killed on hitting 0, that is, 

E x [h{X t ),t <T ] = h(x), t>0,x^0. 

Therefore, we may define a family of measures P* by 

1 



(4.2) 



P£(A) 



■E x [h(X t )l A ,t<T ], x^0,Ae& t , 



for any t 5 s 0. 
(ii) The function h can be represented as 

Pa, (To > e,) 



h(x) 



lim 

giO n(Q > e q 



x # 0, 



where e q is an independent exponentially distributed random variable with para- 
meter q. Furthermore, for any stopping time T and A e and any x 0, 

\imP x (A 1 T<e q \T >e g ) = P|(A). 

This justifies the name 'the stable process conditioned to avoid zero' for the canonical 
process associated with the measures (P|)a.^o- We will denote this process by XK 

Our aim in this section is to prove the following variation of Proposition 14. lT ii) . making 
use of our expression for the density of To, in the stable case only. Our presentation here 
owes much to Yano et al. 33|, §4.3]. 
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Proposition 4.2. Let X be a stable process adapted to the filtration (=^t) t>0 , and h: 



as in 



(i) Define the function 

h(x)n(( > s) 

Then, for any x ^ 0, 

lim Y(s,x) = 1, (4.3) 

s— >00 

and furthermore, Y is bounded on its whole domain. 

(ii) For any x ^ 0, any P x -a.s. finite stopping time T, and A e &t, 

P*(A) = lim P :r (A|To>T + S ). 



Proof. We begin by proving 



= lim P *F°> f , (4.4) 



for x > 0, noting that when x < 0, we may deduce the same limit by duality. 

Let us denote the density of the measure P X (T e •) by p(x, ■). A straightforward 
application of scaling shows that 

P x (T >t) = P 1 (T > x- a t), x>0,t^0, 

and so we may focus our attention on p(l,t), which is the quantity given as p(t) in 
Theorem 13.151 In particular, we have 

p(M) = _^ 2 (Va) sinfrap) T(l - l/a)^ + 0{t - 1/a -^ 
7rsin(7rp) sin(7ra) T(a — 1) 

Denote the coefficient of t 1 /"- 2 j n the first term of this expression by P. Note that all 
limits in Theorem 3.15 exist uniformly in i e [e, oo), so we can integrate in t on any 
interval [to, oo) for t > 0. This allows us to compute P X (T > t) in the asymptotic sense 
from the expansions in Theorem 13.151 by integrating term-by-term. 

To obtain an expression for n(( > t), we turn to Fitzsimmons and Getoor 16J, in 
which the authors compute explicitly the density of n(( e •) for a stable process; see 
p. 84 in that work, where n is denoted P* and ( is denoted R. The authors work with 
a different normalisation of the stable process; they have c = 1. In our context, their 
result says 

/ > , \ a — 1 sin(7r/a) , / 9 , „ . . 

n(C e dt) = —, — ; — r , , ' , ,J ' dt, t ^ 0. 4.5 

vs ; r(l/a) cos(?r(p - 1/2)) V ; 

Denote the coefficient in the above power law by W. 
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We can now compute h. We will use elementary properties of trigonometric functions 
and the reflection identity for the gamma function. For x > 0, 

PsCZo >t) P 



-a-l 



7l(C > t) 



w 



-X 



+ 0{t 



-X 



cos(7r(p — 1/2)) sin(7ra;p) 
T(a) sin(7rp) sin(7ra 
1 sin^pK 
sm(7ra) 

sin(7rap) x 



.a-l 



r(a) 
-rn- 



a 



7T 



+ 0(1). 



This proves (4.4) for x > 0, and it is simple to deduce via duality that this limit holds 
for x ^ 0. 

We now turn our attention to the slightly more delicate result about Y . It is clear 



that the limit in (4.3) holds, so we only need to prove that Y is bounded. We begin by 



noting that, for fixed x 0, Y(t, x) is bounded in t since the function is continuous and 
converges to 1 as t — > go. Now, due to the expression (4.5) and the scaling property of 
X, we have the relation Y(t, x) = Y(\x\ a ~ l t, sgnx). This then shows that Y is bounded 
as a function of two variables. 

With this in hand, we move on to the calculation of the limiting measure. This 
proceeds along familiar lines, using the strong Markov property: 

'P X (1 A ,T > T + s\3? T )' 



P x (A\T >T + s) = E x 



E, 



Ex 



P x (T >T + s) 

Px T (T > s] 



1 ^o>T} Fx{Tq>t + s 

1a1{t >t}KXt)Y(s,X t ) 



n{( > s) 



1 



n(C>T + s) h(x)Y(s + T,x) 



Now, as h is invariant for the stable process killed at zero, (4.2) also holds at T, and 



in particular the random variable h(Xx) is integrable; meanwhile, Y is bounded, as are 
the final two ratios. If we now apply the dominated convergence theorem, we obtain the 
result we seek. □ 

We offer a brief comparison to conditioning a Levy process to stay positive. In this 



case, Chaumont |9|, Remark 1] observes that the analogue of Proposition 4.2(h) holds 
under Spitzer's condition, and in particular for a stable process. However, it appears 
that in general, a key role is played by the exponential random variable analogous to 



that appearing in Proposition 4.1(h) 



4.2 The radial part of the stable process conditioned to avoid zero 

For this section, consider X to be the symmetric stable process, as in section [2j There 
we computed the Lamperti transform £ of the pssMp 

R t = |X t |l{ t<To }, t ^ 0, 
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and gave its characteristic exponent \& in (2.1 
Consider now the process 



R} = \X}\, t>0. 

This is also a pssMp, and we may consider its Lamperti transform, which we will denote 
by £*. The characteristics of the Lamperti-Kiu representation of X$ have been computed 



explicitly in ll|, and the characteristic exponent, of could be computed from this 



information; however, the harmonic transform in Proposition 4.1 (i) gives us the following 
straightforward relationship between characteristic exponents: 

tft(0) = #(0-i( a -l)). 

This allows us to calculate 

r(l/2 - ig/2) r(a/2 + ig/2) 

* w = 2 r((i - a)/2 - ip) me/2) ' ^ EM - 

It is immediately apparent that is the dual Levy process to £. It then follows that 
R is a time-reversal of R*, in the sense of 0, §2]: roughly speaking, if one fixes x > 0, 
starts the process at zero (as in |5], say) and runs it backward from its last moment 
below some level y, where y > x, simultaneously conditioning on the position of the left 
limit at this time taking value x, then one obtains the law of R under P x . 

We remark that this relationship is already known for Brownian motion, where R} is 
a Bessel process of dimension 3. However, it seems unlikely that any such time- reversal 
property will hold for a general Levy process conditioned to avoid zero. 

4.3 The entrance law of the excursion measure 



It is known, from a more general result [13J, (2.8)] on Markov processes in weak duality, 
that for any Borel function /, 



I™ e-« t n(f(X t ))dt= f f{x)E x [e-" To ]dx. 

JO JR 



(This formulation is from 33|, (3.9)].) In terms of densities, this may be written 



n(X t e dx) dt = P X (T e dt) dx 

= |x| _a P sgn( _ :E )(ro e \x\~ a dt)dx 

Therefore, our expressions in Theorem 13.151 for the density of T yield expressions for 
the density of the entrance law of the excursions of the stable process from zero. 

4.4 Identities in law using the exponential functional 



In a series of papers (Bertoin and Yor [3|, Caballero and Chaumont [5(, Chaumont et al. 
0) it is proCed that under certain conditions, the laws (PJ„„ of an a-pssMp X admit 
a weak limit Pq as x J, 0, in the Skorokhod space of cadlag paths. If £ is the Lamperti 
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transform of X under Pi, then provided that E|£i| < go and m := E£i > 0, it is known 
that the entrance law of Po satisfies 

E (/(Xf)) = -J-Eili-aO-'fit/H-aO)), 
am 

for any t > and Borel function /. Similar expressions are available under less restrictive 
conditions on £. 

It is tempting to speculate that any rssMp may admit a weak limit Po along similar 
lines, but we do not propose any results in this direction; instead, we demonstrate similar 
formulae for the entrance law n{X t s •) of the stable process, and the corresponding 
measure Pq for the stable process conditioned to avoid zero. 

Let X be a stable process, possibly asymmetric. From the previous subsection, we have 
that 

n(f(X t ))= \x\~ a p(sgn(-x), \x\~ a t)f(x) dx. 

J— 00 

Substituting in the integral, and recalling that the law of T Q for the stable process is 
equal to the law of the exponential functional I(a£) of the Markov additive process 
associated with it, we obtain 

"00 



n(f(X t )) = - I p(l,«)/(-(t*A)- 1/ >- 1/ot t 1/o, - 1 d« 

r>oo 



a 



a 



+ - p(-l,u)/((u/t)" 1/a )«" 1/ai t 1/ai " 1 clti 



^[/(-(t/JM^/Kr 1 ^- 1 ] 

a L J 

+ -E 2 [f((t/I(aO) 1/a )I(aO- 1/a t 1/a - 1 ]. 
a 



Recall from 27] that the law Pq of the stable process conditioned to avoid zero is 
given by the following harmonic transform of the stable excursion measure n: 

Pj(A) = n(l A h(X t ),t < C), t > 0, A e & u 



with h as in (4.1) Therefore, applying the above result to the Borel function hf, we 
obtain 

El(f(X t )) = n(h(X t )f(X t )) 

= r(-a)^^E 1 [j(<)- 1 /(-(^K)) 1/a )] 

7T 

+ r(-a)^^E 2 [j(<)-V((V/(«C)) 1/a )], 

where we emphasise that I(a£) (under Ej) is the exponential functional of the Markov 
additive process associated to X. 
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